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The effect of a number of mechanisms designed to suppress decoherence in open quantum systems 
are studied with respect to their effectiveness at slowing down the loss of entanglement. The 
effect of photonic band-gap materials and frequency modulation of the system-bath coupling are 
along expected lines in this regard. However, other control schemes, like resonance fluorescence, 
achieve quite the contrary: increasing the strength of the control kills entanglement off faster. The 
effect of dynamic decoupling schemes on two qualitatively different system-bath interactions are 
studied in depth. Dynamic decoupling control has the expected effect of slowing down the decay 
of entanglement in a two-qubit system coupled to a harmonic oscillator bath under non-demolition 
interaction. However, non-trivial phenomena are observed when a Josephson charge qubit, strongly 
coupled to a random telegraph noise bath, is subject to decoupling pulses. The most striking of 
these reflects the resonance fluorescence scenario in that an increase in the pulse strength decreases 
decoherence but also speeds up the sudden death of entanglement. This demonstrates that the 
behaviour of decoherence and entanglement in time can be qualitatively different in the stong- 
coupling non-Markovian regime. 

PACS numbers: 03.65.Yz, 03.67.Bg, 03.67.Pp 



I. INTRODUCTION 



Entanglement is one of the basic features that distin- 
guish quantum systems from their classical counterparts, 
and has its origins in the inherent non-locality of quan- 
tum mechanics [l| . It is the most useful resource in quan- 
tum information theory and is indispensible for di- 
verse quantum information tasks such as quantum com- 
munication, teleportation, quantum error correction, su- 
perdense coding, one-way communication etc. In closed 
systems — that is, systems which are completely isolated 
from their surroundings — entanglement remains con- 
served under a local unitary evolution, and decays slowly 
under a nonlocal evolution. This makes these systems 
ideal for quantum information tasks. Closed systems are, 
however, a rarity in the natural world. More often than 
not, quantum systems are open, that is, they are in con- 
tact with the surrounding environment — a thermody- 
namic reservoir, for example Quantum systems 
are extremely fragile, and the dissipative effects of the 
environment gives rise to the phenomenon of quantum 
decoherence Q. As a result, the system undergoes an 
asymptotic transition to classicality and hence loses all its 
entanglement, which is a purely quantum phenomenon. 
This happens even if the evolution is a local unitary one. 
Nevertheless, this in itself is not a bad scenario, for if the 
decoherence rate is low, then entanglement takes a long 
time to completely disappear and such systems can func- 
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tion as useful quantum devices for sufficient periods of 
time. However, recent studies 0, II] have uncovered sys- 
tems where the rate of loss of entanglement is exponen- 
tially higher than the decoherence rate. This results in a 
finite time to classicality, and consequently, a finite time 
to the total loss of entanglement — a phenomenon given 
the name entanglement sudden death (ESD). Systems 
that suffer from ESD are rendered unusable for quan- 
tum tasks. Naturally then, ESD has dire implications 
for the success of quantum tasks, and has become one 
of the premier branches of quantum information study 
in recent times. Some of us have recently investigated 
this phenomenon for the case of n-qubit states at finite 
temperature [§| , as well as for spatially separated n-mode 
Gaussian states coupled to local squeezed thermal baths 

Given the obvious importance of ESD regarding the 
success of quantum tasks, it is thus a worthwhile exer- 
cise to investigate ways and means of controlling the rate 
of loss of entanglement. Error-correcting codes [ll|, [i3 
and error-avoiding codes [l3| ( which are also knows as 
decoherence-free subspaces) are such attempts. Open 
loop decoherence control strategies [l^^Tsj are another 
class of widely used strategies used to this effect, where 
the system of interest is subjected to external, suitably 
designed, time-dependent drivings that are independent 
of the system dynamics. The aim is to cause an effective 
dynamic decoupling of the system from the ambient envi- 
ronment. A comparative analysis of some of these meth- 
ods has been made in [l9| . Another mechanism known to 
slow down the process of decoherence is through manip- 
ulation of the density of states. This has been put to use 
in photonic band-gap materials, which is used to address 
questions related to the phenomenon of localization of 
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light [13421. 

In this paper, we analyze the evolution of entangle- 
ment in two-qubit systems connected to local baths (or 
reservoirs). A number of studies of entanglement in open 
quantum systems have been made [25,-,27|. Here we ad- 
dress the need to have a control on the resulting nonuni- 
tary evolution, as motivated by the above discussion, and 
study several methods of doing so. These include ma- 
nipulation of the density of states in photonic crystals, 
modulation of the frequency of the system-bath coupling 
and modulation of external driving on two-qubit systems 
as examples of systems undergoing Markovian evolution. 
We also study control methods in systems undergoing 
non-Markovian evolution. The first of these is dynamic 
decoupling — which is an open-loop strategy — on a 
two-qubit system that is in contact with a harmonic os- 
cillator bath. This system undergoes a quantum non- 
demolition interaction, where dephasing occurs without 
the system getting damped. The second is a Joscphson- 
j unction charge qubit subject to random telegraph (1//) 
noise due to charge impurities. The plan of the paper 
is as follows. In Section II, we introduce the basic tech- 
niques and formalism used in this paper, including the 
formal way of solving the Lindblad master equation. We 
also introduce the phenomenon of channel-state duality 
and the law of entanglement decay, both of which will be 
used subsequently. In Section III, we study the evolution 
of entanglement in photonic band gap materials and the 
effect of the special characteristics of such materials on 
ESD. The effect of frequency modulation of the system- 
bath coupling on ESD is studied in Section IV. This is 
followed by a study of the effect of ESD on a two-level 
system driven by an external resonant field which is in 
resonance with the transition frequency. Finally, in sec- 
tion VI, we study the effect of dynamic decoupling on the 
evolution of entanglement and ESD. We pay particular 
attention to the Bang-Bang strategy [H with regard to 
the usual two-qubit system under a QND interaction in 
section VI (A) and also to a Josephson-junction charge 
qubit subjected to random telegraph noise, and make 
comparisons. We conclude our paper in section VII with 
further discussions. 



II. PRELIMINARIES 

An open quantum system, as defined in the introduc- 
tion, is exposed to its environment, which is usually a 
thermal reservoir. The dynamics of such a system is nat- 
urally dictated by its interaction with its environment. 
If H be the total Hamiltonian of an open system, then 
H = Hs + Hb -I- HsR, where Hs and Hn are the system 
and reservoir Hamiltonians respectively and Hsr is the 
interaction Hamiltonian. Open systems undergo nonuni- 
tary evolution due to this interaction term, and, depend- 
ing on the type of the system-reservoir {SR) interaction, 
can be broadly divided into two categories — dissipativc 
and non-dissipative. In the former, the system Hamilto- 



nian does not commute with the interaction Hamiltonian, 
[Hsi Hsr] 0, and dephasing occurs along with dissipa- 
tion and decoherence. In the latter however, these two 
do commute — [Hs,Hsr] =0 — and hence the SR in- 
teraction is characterized by a class of energy-preserving 
measurements where dephasing occurs without damping 
the system [2^ [2^. Such a non-dissipative system, as 
well as the corresponding interaction, is called a Quan- 
tum Non-Demolition (QND) system. 

We are interested in the time evolution of OQS, i.e., of 
the system state ps- Let the initial state of the system- 
bath combination be p(0), and let the state at time t be 
p{t) = Up{0)U\ where U = e~''^* is the time evolution 
operator. The state of the system alone is obtained from 
p{t) by simply tracing out the bath degrees of freedom: 
psit) — Tr^j [p(t)]; where Ttr implies a partial trace over 
the bath. The evolution of the system-bath combination 
in unitary, and 

m = -^[H,pit)] (1) 

is the equation of motion. However, the evolution of 
the system itself is nonunitary, and thus requires a more 
general equation of motion which, after the application 
of the Born, Markov and rotating wave approximations, 
can be written as 

ps{t) = ^t[Hs,ps{t)]+ (2) 

Y.^, (f.PsWf] - ^-{ps{t),FjF,}y 

This is a master equation in its Lindblad form. It can be 
written in super operator form as 

ps{t) ^ C[ps{t)] (3) 

Psijit) ^^^^^L^j^kiPskiit), (4) 
fe, / 

where C is the super operator acting on the system state 
Psit) and is effectively a time-derivative, and where L 
is the matrix representation of £. In general, C is time 
independent and the solution of the above equation can 
be written formally as 

ps{t) ^ A[ps{0)] (5) 
Psijit) = YVijMPski{0), (6) 

k,l 

where V = exp(Li) is the matrix representation of the 
time evolution map A. If the system is evolving under 
unitary evolution w then the matrix V is simply w<Siw*, 
where w* represents the complex conjugate of w in a fixed 
basis. 

Channel-State Duality- A quantum channel is a con- 
duit for the transmission of quantum as well as classical 
information, and is essentially a completely positive map 
between spaces of operators. Any such physical quan- 
tum channel acting on a d-dimensional quantum state 
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can be mapped to a positive operator in (f dimensions, 
and, if the channel is trace-preserving, then the corre- 
sponding positive operator will have unit trace. Similarly, 
a valid density matrix in d^-dimensions can be mapped 
to a trace preserving physical channel in d dimensions. 
Such a two-way mapping between a quantum state and 
its corresponding higher dimensional operator is called 
channel-state duality. 

The time evolution operator A is a physical quantum 
channel represented by the matrix V. If M be a valid 
density matrix corresponding to the map A, it is given 

by mM\ 



M = {I(g)A)[\(t)+){<p^ 



(7) 



where \(f>'^} = X^iLi N) ^ maximally entangled state 
in d (E) d. Here the channel is applied to one side of a 
maximally entangled state in d ® d. We shall use the 
symbols V and M to represent the matrix representation 
of the time evolution map A and a valid density matrix 
corresponding to it, respectively, throughout the paper. 

Law of entanglement decay This law says that the 
evolution of entanglement in a bipartite entangled state 
under a local one-sided channel can be fully character- 
ized by its action on a maximally entangled state. The 
amount of entanglement at any time t in a given initially 
entangled two-qubit state |x), under the action of a one- 
sided quantum channel, is equal to the product of the 
initial entanglement in the given state and the entangle- 
ment in the state which we get by applying the channel 
on one side of a two-qubit maximally entangled state. 
Mathematically this can be written as: 

C((I®A)[|x)(xl])=C([|x)(x|])C((I®A)[|0+)(0+|]), 

(8) 

C(.) being the concurrence [111. Therefore, it is enough 
to study the evolution of state. 



III. EVOLUTION OF ENTANGLEMENT IN 
THE PRESENCE OF PHOTONIC CRYSTALS 

In this section wc consider a system of two level atoms 
interacting with a periodic dielectric crystal. The peri- 
odic dielectric structure gives rise to the photonic band 
gap [l^, m, [sJl . The effect of this particular structure 
on electromagnetic waves is analogous to the effect semi- 
conductor crystals have on the propagation of electrons, 
and leads to interesting phenomena like strong localiza- 
tion of light [131 , inhibition of spontaneous emission 
and atom-photon bound states [l^, [H, . The origin 
of such phenomena can ultimately be traced to the pho- 
ton density of states changing at a rate comparable to 
the spontaneous emission rates. The photon density of 
states are of course estimated from the local photon mode 
density which constitutes the reservoir. 

The advantage of studying them is that the presence of 
the photonic gap suppresses decoherence [s^l . Let us con- 
sider a two-qubit system. One of these is locally coupled 



to a photonic crystal reservoir initially kept at tempera- 
ture T = 0. In this case the entanglement dynamics can 
be obtained by studying any one of the qubits individu- 
ally. We start with the following Hamiltonian: 

H = ycr^ + ^ujkalak + ^{gkal(T^ + glaka+), (9) 
fe fc 

where wq is the natural frequency of the two level atom, 
Ldk is the energy of the kth mode and gk is the frequency 
dependent coupling between the qubit and the photonic 
crystal, which is acting as the reservoir here. If we restrict 
the total atom-reservoir system to the case of a single 
excitation , the evolution of a given state of the qubit 
is then given by [STf : 

f Poo{OMt)\' Poi(0)c(t) 



(10) 



where 

c{t) =e f A+e^^+*(l + <^>{X+e"/^^m 



-A- 



iX~_t 



(1 + $(A_ 



in/i 



$(a;) = 



OO 



,ii5t 



is the error function 



_4j' 



AS 



2 



^0 



Here (5 = wq — Wc is the detuning of the atomic frequency 
and ujc is the upper band-edge frequency. We have made 
use of the following photon-dispersion relation near the 
band edge: « Wc + A{k — fco)^, where A « ujc/kQ, d is 
the atomic dipole moment and eq is the vacuum dielectric 
constant. 

The density matrix p{t) in (fTU]) is related to the initial 
density matrix p(0) by a map A, given by p{t) — A[p(0)], 
whose matrix representation is 

/ 

c(t) 


Vi-|c(t)P 



V 




(11) 



Channel-state duality, explained earlier in Section II, en- 
sures that there exists a two-qubit density matrix M for 
every single-qubit channel V. This matrix M can be 
written as 



M = Ay (!$+)($+ 1) 






c{t) 



\c{t)\' 

1 



(12) 
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where |<I>+) = (|00) + |11)) is a two-qubit maximally 

entangled state. The concurrence of M is |c(<)p. The 
complex function c is a function of the detuning param- 
eter 6 and time t. Therefore, we need to see the effect of 
a on entanglement. If we assume that 5 = Aa^, c can 
then be written in the following simplified form; 



c{t) 



VI -4A2 



i(d+e<^(l + $(d+e-/4V7)) 



where d± = — 1 ± \/l — 4 A and r = a^t. Since the 
entanglement in M is |c(t)p, it is now a function of 6 
and T. 

The evolution of entanglement for different values of 
A can be seen in the FIGS. [I(a)| [1(b)] The inset of the 
figures depict the evolution of entanglement, computed 
using concurrence (see appendix), for the usual case of 
zero band gap, while in FIGS. 1(a) l(b)[ evolution of en- 
tanglement is shown for increasing influence of the band 
gap. In FIG. l(a)[ the system is within a gap in the pho- 
tonic spectrum, indicated by the negative value of A and 
hence also (5, as a result of which coherence is preserved 
and the decay of entanglement is arrested. This feature 
is further highlighted in FIG. 1(b) which is also for the 



case of negative A of higher order of magnitude than that 
in FIG. 1(a) and as a result there is a greater persistence 
of entanglement. Thus we find that with the increase in 
the infiuence of the photonic band gap on the evolution, 
entanglement evolution is preserved longer. From Eq. 
([T^. it can be seen that, following the arguments of the 
previous section, there is no ESD in this feature 
corroborated by the FIGS. [TJ 



IV. FREQUENCY MODULATION 

Agarwal and coworkers (isj introduced an open- 
loop control strategy which involved modulation of the 
system-bath coupling, with the proviso that the fre- 
quency modulation should be carried out at a time scale 



which is faster than the correlation time scale of the heat 
bath. The technique of frequency modulation has been 
used earlier to demonstrate the existence of population 
trapping states in a two- level system (39| . Raghavan et 
al. [40[ showed the connection between trapping in a two- 
level system with frequency-modulated fields in quantum 
optics and dynamic localization of charges moving in a 
crystal under the action of a time-periodic electric field. 

Consider the Hamiltonian given in Eq. ([9]) . Frequency 
modulation essentially involves a modification of the cou- 
pling g — the modulated coupling is gexp{— MTisinz^i}, 
where m is the amplitude and v is the frequency of the 
modulation. Choosing m such that Jo{m) = 0, where Jq 
are the Besscl functions of order zero, the decay of the 
excited state population can be shown to be significantly 
arrested. The resulting master equation in the interac- 
tion picture, when a ppl ied to the evolution of a standard 
two- level system, is [18l |: 



dp 



{C^+iS+S-p- S-pS+) 



dt [k - iA)2 + zy2 

+C+-{S-S+p - S+pS~)] + h.c, A^{uo- cj). 

(13) 

where the Besscl function expansion g^-tn^smiut) _ 
Y^\=°°oo Ji{'m)e~^^^'^ has been used. Additionally, the 
modified bath correlation functions are assumed to have 
the forms C-+{t) = Co~+e-''*e'"* and C+-{-t) = 
where k is the bath correlation frequency. 
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Now, we have 
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=>p(t) ==exp(At)p(0), (15) 
=^ exp(At),j,fcip(0)fc( = V{t)tjMp{^)ki, 

(16) 
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where V = exp(Lt) and L is the matrix representation 
of A. We obtain the matrices L and V using Eq. ([Tc 
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/ i (Co-+e-2R°(")^* + C+-) 



V = cxp(Lt) 



-2Re(a)C^- / 
^(Co+-(l-e-2R°(«m)) \ 
e-"'^* 

e-"'" 

y ^(C-+(l-e-2R°("m)) i (Co+-e-2R°(«m + Co-+) y 

I 



(17) 
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where a = %Z\^y\^^^ and T = + C+- . If A/ 



(I® X^)(|0+)(0+|), then we have 



M 



( Mil e-"^* \ 

M22 

M33 

\e-"*" M44 J 



(19) 
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where 



Mil 
M33 

Af44 



[I _ g-2Ro(a)Tt^| 
(1 _ g-2Rc(a)Tt^ 



(20) 
(21) 
(22) 
(23) 



If this state M is separable, the factorization law for en- 
tanglement decay [S^] allows us to assert that all states 
will show ESD. The state M will be separable only if the 
condition 



2X 



J.2 



(24) 



holds, where X = exp(— 2Re(a)r<:). The roots of the 
above equation are 




(25) 



The negative root is less than unity (X_ < 1), implying 
that there exists, always, a finite and positive time tESD 
at which the system loses all its entanglement. This is 
given by 



tESD 



1 



2Rc(a)T 



log(X. 



(26) 



The modulation factor v appears in the numerator and 
therefore it can be expected that a higher frequency of 
modulation should sustain entanglement longer. This is 
confirmed in the plot of Iesd against v (FIG. [2]). This 
result is not altogether surprising, for a higher degree of 
modulation is naturally expected to filter out the influ- 
ence of the bath and increase the coherence which ul- 
timately results in entanglement sustaining for a longer 
period of time. 



V. RESONANCE FLUORESCENCE 

In the previous section, wc focused on the decrease in 
the time to ESD through an increase in the degree of fre- 
quency modulation of the system-bath coupling. In this 
section, we study a system where a two-level atomic tran- 
sition is driven by an external coherent single-mode held 
which is in resonance with the transition itself. We shall 
show that, in this situation, an increase in the Rabi fre- 
quency — which plays the role of the modulator — pro- 
duces the opposite effect by speeding up ESD. The behav- 
ior of such driven systems has been well studied in the lit- 
erature and has found many applications. In contrast to 
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(a) Entanglement (concurrence) of M as a function of r 
for detuning parameter A = —0.1. In this plot the 
behaviour of entanglement is different from the one in the 
case for A = (see inset). Here entanglement is seen to 
converge to a non-zero value at large r. 
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(b)Entanglement (concurrence) of M as a function of r 
for detuning parameter A = —0.25. This plot shows that 
higher the detuning , larger will be the asymptotic value 
of entanglement. 

FIG. 1: Evolution of entanglement in photonic band gap crys- 
tals at zero temperature. 



the situation here. Lam and Savage [4l| have investigated 
a two-level atom driven by polychromatic light. The phe- 
nomenon of tunneling in a symmetric double- well poten- 
tial perturbed by a monochromatic driving force was an- 
alyzed by Grossmann et al., (42| . while photon-assisted 
tunneling in a strongly driven double-barrier tunneling 
diode has been studied by Wagner [i^ . 

The analysis of the said driven system begins with its 
Hamiltonian which, when written in the interaction pic- 
ture, is HsR = -E{t)-D{t). RcTcE{t) = ee-^"°*-he*e*'^«* 
is the electric held strength of the driving mode (treated 
classically), luq is the atomic transition frequency and 
D(t) is the dipole moment operator in the interaction 
picture. The driven two-level system is coupled to a 
thermal reservoir of radiation modes. If 70 be the spon- 
taneous rate due to coupling with the thermal reservoir 
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FIG. 2: In this plot we have time of BSD, i.e, the time at which 
a maximally entangled initial state loses all its entanglement 
and become separable when exposed to a bath, against the 
frequency of modulation v. In this case we have kept the 
value of m to be the first zero of the Bessel function Jo, i.e, 
m = 2.4048. Also, at i/ = the value of Iesd is 1.4. The 
value of the other parameters are: k — 0.1, A = 0.1, Cq~ = 



0.1. 



and N = A^(a;o) be the Planck distribution at the atomic 
transition frequency wq, then the evolution of this com- 
posite system is given by the following master equation 



7o(A^ + l 



2 

70 (TV) 



[2(T_p(t)(T-|_ — a^a^p[t) — p{t)a+a- 

[2cr+p(t)CT_ - a^a+p{t) - p(t)cr_cr+] , 

(27) 



where Q. — 2e -d* is the Rabi frequency, d being the tran- 
sition matrix element of the dipole operator. The term 
— (ri/2) [cr+ + <T_] characterizes the interaction between 
the atom and the external driving field in the rotating 
wave approximation. As usual, <7± are the atomic rais- 
ing and lowering operators, respectively. 

Let us take two identical qubits and, as before, assume 
that one of them interacts locally with a thermal bath 
and is subject to monochromatic driving by an exter- 
nal coherent field. The master equation (Eq. Wf) yields 
the corresponding matrices Vrf and M^f (where the sub- 
script rf stands for resonance fluorescence): 



Vrf = 



Mrf = -: 




(28) 



(29) 



ai + a4 =1 + - X'^ \ cos{pt) - ^ sin(/it) ) ) 6*3 
+ sin(^i) + 
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cos(/it) — — sin(/it) 



02 = — X^smlypit), 



in 



61 + 64 = - X'\S+ + S-)- —X^sm{pt)S3 

+ X^ (^cos{pt) + J- sin(^i)^ {S^ - S+) 
ifl 

61 - 64 = — X^sin(yut), 
/i 

&2,3 ^^X^ ± X^ (^cos{pt) + J- sm{pt) 



di + d4 ~2— (fli + 04), 
di — di ~ — (ai — 04), 



It 

X =e--. 



»^7o 

72 + 2n^ 



7o7 



72 -I- 21]2 ^ 
7=7o(2A^+l), 



/i =Vn^ (7/4)2. 

Using these, we plot, in FIG. ([3]), concurrence vs the time 
to ESD for different values of the Rabi frequency il and 
observe that tssD decreases for an increase in fl. This 
is contrary to the result derived in the previous section, 
where an increase in the modulation frequency v delayed 
the loss of entanglement. The decrease in tESD does not 
however continue indefinitely, but rather saturates to a 
certain value for large values of the Rabi frequency. This 
is an interesting result, further analysis of which will be 
carried out in a future work. 

Let us now consider the situation where the system, 
consisting of the excited two-level atom, is at zero tem- 
perature. Let us also consider the evolution of entangle- 
ment for two cases demarcated by the relation between 
the Rabi frequency and the spontaneous rate of coupling 
with the thermal reservoir. For the underdamped case 
when n > 7o/4, the quantity p is real (since = 
at T = 0) and hence both the upper level occupation 
and coherence exhibit exponentially damped oscillations. 
Conversely, in the ovcrdamped case, ft < ^ ^ p G I 
and both these quantities decay monotonically to their 
stationary values. The evolution of entanglement, how- 
ever, works in an opposite way . Entanglement decays 
faster for the underdamped case than for ovcrdamping, 
where the Iesd is higher. One possible reason for this 
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FIG. 3: Entanglement (concurrence) of Mrf as a function of 
time for different values of Rabi frequency Q, varying from 
to 0.5: (the last curve on the right hand side) corresponding 
to pure damping and 0.5 (first curve on the left hand side) 
corresponding to the underdamped case, i.e., it covers both 
the overdamped as well as the underdamped cases. In the 
inset one can see that as we increase the fl the tssD seem to 
converge a,t t = 17.0. Here 7 = 0.1. 



could be relative positions of the three Lorentzian peaks 
of the inelastic part of the resonance fluorescence spec- 
trum. The central peak is at w = luq and the rest are 
at w = cjQ ± /i for the underdamped case, whereas all 
three peaks are at a; = cjo for the overdamped case. This 
indicates that the decay of entanglement in the under- 
damped should be closely dependent on the quantity /i. 
This in turn depends on both the dissipation parameter 
7 and the Rabi frequency, the latter in itself a function of 
the driving strength of the external field and the dipole 
transition matrix elements. Thus, in the underdamped 
case, there exists greater avenues for the decay of quan- 
tum coherences as well as entanglement than the over- 
damped case. Phenomenologically, for the underdamped 
case (f2 > 7o/4), the two-level atom interacts with the 
external monochromatic field multiple times before spon- 
taneously radiating a photon (see, for example, chapter 
10 of ref. nil)- Such numerous interactions allows quan- 
tum correlations to develop between the two atomic levels 
and the quantized levels of the field. The phenomenon 
of monogamy of entanglement (45| thus ensures that the 
amount of quantum correlation between the two qubits 
will decrease. Additionally, it can be seen that at a higher 
Rabi frequency, 51 dominates the dissipation and thus 
causes a saturation of the time to ESD, as shown in FIG. 

m 



VI. DYNAMIC DECOUPLING AND THE 
EFFECT ON ESD 

As discussed in the Introduction, open-loop control 
strategies involve application of suitably tailored control 



fields on the system of interest, with the aim of achieving 
dynamic decoupling of the system from the environment 
Bang-Bang control is a particular form of 
such decoupling where the decoupling interactions are 
switched on and off at a rate faster than that set by the 
environment. The application of suitable radio frequency 
(RF) pulses, applied fast enough, averages out unwanted 
effects of the environment and suppresses decoherence. 
In this section, we compare the effect of Bang-Bang de- 
coupling on the evolution of entanglement in systems con- 
nected to two different types of baths. One bath type 
is composed of infinitely many harmonic oscillators at a 
finite non-zero temperature T and couples locally to a 
two- level atom acting as the qubit, while the other adds 
random telegraph noise to a Josephson-junction charge 
qubit. It has been shown for the former case that all 
two-qubit states shows ESD at finite T 



A. Bang-Bang decoupling when the bath consists 
of harmonic oscillators 

1. Quantum Non-Demolition Interaction 

Let us consider the interaction of a qubit with a bath of 
harmonic oscillators where the system Hamiltonian com- 
mutes with the interaction Hamiltonian so that there is 
no exchange of energy between the system and the bath 
— this is quantum non-demolition dynamics [28l . [29j . 
The only effect of the bath will be on the coherence el- 
ements of the qubit evolution, which will decay in time 
at the rate 7. The total Hamiltonian for the system plus 
bath is: 

Ha = H, + Hb+Hi- (30) 

k 

Hi = ^(^z{gkbl+ glhk). 

k 

Here the system Hamiltonian Hq commutes with the in- 
teraction Hamiltonian Hi and the evolution of such a 
system is called pure dephasing. For simplicity we will 
work in the interaction picture where the density matrix 
of the system plus bath and the interaction Hamiltonian 
transform as: 

p(i) = e*(ff.+«i^)*p(f)e-^(»^+^«)*, (31) 
Hit) = <J, ^(f'^^I^"''* + dkhe-'"^"'). (32) 

k 
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From here we can write the total time evolution operator 
for the system plus bath as 



U{to,t) = Tcxp <j / dsH{s) 

J to 



cxp 



(33) 



where ^k(t) = —(1 — exp(icjfct)). We are interested in 
calculating 

~poi{t) = (OlTrs [u{h,t)pm)U\to,t)^ (34) 

Assuming that in the beginning the bath and the qubit 
were uncorrelated and that the bath is in a thermal state, 
we have [l^ : 



Poi(i) =Poi(to)e-^(*'"*\ 



where 



coth(^) 



(35) 
(36) 



The matrix representation of the evolution operator V 
can be written from here as: 



V 



/I {)\ 

e-7(to,t) 

e-T(*0'*) 

V 1 / 



(37) 



To see the evolution of entanglement, it is enough to 
evolve the maximally entangled state 10"'") = (1 00) + 
|ll))/-\/(2) and study the evolution of entanglement. 
Evolving one subsystem in state | </)"'") gives rise to the 
density matrix: 



M 



( 



V 



1 




-l(to,t) 



e-7(*o,t) 




1 



(38) 



The concurrence in the state M is directly proportional 
to e""^(*«^*). 



2. Dephasing under Bang-Bang dynamics 

The aim of Bang-Bang decoupling, acting on the sys- 
tem of interest, is to hit it with a sequence of fast radio- 
frequency pulses with the aim of slowing down decoher- 
ence (see FIG. S]). Adding the radio frequency term to 
the system-plus-bath Hamiltonian iJo (Eq. ((30)) ). we get 



Tl=l 



FIG. 4: Sketch of the pulse sequence used in bang-bang de- 
coupling procedure. 



where £ 



(n) 
V 



to + nAt, n = 1, 2, • • • , and 



C/(")(i) 



U t^p^ <t<t 
elsewhere. 



(n) 



(41) 



The term Hijp acts only on the system of interest — 
here the qubit. It represents a sequence of Up identical 

pulses, each of duration Tp, applied at instants t = 
The separation between the pulses is r = At. The decay 
rate for this pulsed sequence evolution is |14| : 



jpiN,At)^J2 



\Vk{N,At)\^ 



coth(^) 



(42) 



H{t) = Ho + HRFiL^o,t), (39) 

rip 

HRF{t) = J2 U^"^{cos[ujoit - 4"))]a, + sin[c^o(t - t^"^)]^y 

(40) 



where 

|r;fc(7V,L^feA0P = 4(1 - cos(c^;,At))2 

/ N-l \ 

X iV+ E 2ncos[2{N - n)ujkAt] . 

\ n=0 J 

(43) 

In [iJl it has also been shown that jr/fcp < j^/jp which 
implies that decoherence is suppressed. Also, it is evi- 
dent that a lower value of rj implies a lower value of 7. 
Consequently, we can conclude that that Bang-Bang de- 
coupling slows down entanglement decay. 

B. Josephson Junction qubit 

Although solid state nanodevices satisfy the require- 
ments of large scale integrability and flexibility in de- 
sign, they are subject to various kinds of low-energy ex- 
citations in the environment and suffer from decoherence 
problems. There have been a number of proposals in 
this context about the implementation of quantum com- 
puters using superconducting nanocircuits [i^ \^ . Ex- 
periments highlighting the quantum pro perties of such 
devices have already been performed |48l . |49| . Here the 
concept of a Josephson-junction qubit comes into promi- 
nence. A charge- Josephson qubit is a superconducting is- 
land connected to a circuit via a Josephson junction and a 
capacitor. The computational states are associated with 
charge Q in the island and are mixed by Josephson tun- 
neling. For temperatures much lower than the Josephson 



energy, ksT <^ Ej |50l - l52l |. we have the Hamiltonian 



He 



e 
2 



E, 



(44) 
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with the charging energy Ec dominating the Josephson 
energy. Here, e = e{V) = 4:Ec{l - C2V/e), C2 is the 
capacitance of the capacitor connected to the island and 
V is the external gate voltage. 

An important source of decoherence in the operation of 
Josephson charge qubits is due to fluctuating background 
charges (BCs) (charge impurities) . These are believed to 
originate in random traps for single electrons in dielectric 
materials surrounding the superconducting island. These 
fluctuations cause the 1// noise, also known as random 
telegraph noise, at low frequencies, and are directly ob- 
served in single electron tunneling devices [s^ [s^l • This 
has also been studied in the context of fractional statistics 
in the Quantum Hall Effect ISSll . This noise, arising out 
of decoherence, is modeled jSll . Is^ by considering each 
of the BCs by a localized impurity level connected to a 
fcrmionic band, i.e., the quantum impurity is described 
by the Fano Anderson model. This is the quantum ana- 
logue of the classical model of TV independent, randomly 
activated bistable processes. For a single impurity, the 
total Hamiltonian is: 



(45) 



where 



H' = e.b^b + [Tkclb + h.c^ + J2 ^kclck. (46) 

fe k 

Here describes the BC, b represents the impurity 
charge in the localized level ec, Ck the electron in the 
band with energy e^, and Hq is as in Eq. (|44|) . The im- 
purity electron may tunnel to the band with amplitude 
Tfc. The BC produces an extra bias v for the qubit via 
the coupling term {v/2)b^baz- An important scale is the 
switching rate 7 = 27rp(ec)|Tp, where p(ec) is the den- 
sity of states of the band and it is assumed that we are 
working in the the relaxation regime of the BC where the 
tunneling rate to each fermionic band is approximately 
equal. The fraction v/j determines whether the opera- 
tional regime of the qubit is weak (w/7 ^ 1) or strong 
(v/7 > 1). Studying the single BC case is important, 
since it has been shown 5l| that the effect of multiple 
BCs can be trivially extended from that of a single BC. 
For multiple strongly coupled BCs producing 1// noise, 
the effect of a large number of slow fiuctuators is minimal 
and pronounced features of discrete dynamics such as sat- 
uration and transient behavior are seen. There are two 
special operational points for the qubit (|44l): (a) e = 0, 
corresponding to charge de gen eracy and (b) Ej = 0, for 
the case of pure dephasing |56l . [STj , where tunneling can 
be neglected. We will consider this case later in detail 
and make a comparison of ESD, for the case of pure de- 
phasing, between the harmonic oscillator and 1// baths. 

In order to study the effect of the BC on the dynamics 
of the qubit, the general procedure is to calculate the uni- 
tary evolution of the entire system plus bath and trace 
out the bath degree of freedom, i.e, pgit) — tTE{W{t)}, 
W{t) being the the full density matrix. In the weak cou- 
pling limit a master equation for PQ{t) can be written 



[58l | . The results in the standard weak coupling approach 
are obtained at lowest order in the coupling w, but it has 
been pointed out that higher orders are important for a 
1// noise 13,1511. 

The failure of the standard weak coupling approach is 
due to the fact that the 1// environment includes fiuctu- 
ators which are very slow on the time scale of the reduced 
dynamics. To circumvent this problem one considers an- 
other approach in which a part of the bath is treated on 
the same footing as the system [s^l- We study the evo- 
lution of this new system and later trace out the extra 
part which belongs to the bath, i.e., p{t) = Trfb{PF(i)}. 
We then obtain pQ{t) from p{t) as PQ{t) = Trb{p(i)}, 
where fb stands for fermionic band. In that regard we 
split the Hamiltonian (|^5|) into a system Hamiltonian 

i/o = E[q — — b'^baz + ecb^b and environment Hamiltonian 

= J2k '^kc\ck coupled by V = X]fe Tkc\b -|- hj 

The eigenstates of Hq are product states of the form 
|0)|n),e.g, 

|a) = |0+)|O), 

\b) = \e-m, 

with corresponding energies 

"2"' 2"' ^^^^"^ y^^"- 

Here 16*-!-) are the two eigenstates of (Jn, the direc- 
tion being specified by the polar angle 6 and (p = 

0. The two level splittings are f2 
n' = ^{e + v)^ + E^, and cos(6l) 



E^ and 



J, sin(6l) = 

E^/n, cos(0') = {e + v)/n', sin(6l') = Ej/Q'. 

In the basis of the eigenstates of Hq the master equa- 
tion for the reduced density matrix p{t) in Schrodinger 
representation reads: 



dpijjt) 
dt 



= -iuj,pij{t) + 



mn 



R 



it), 



(47) 



where uJij is the difference of the energies and Rij^mn are 
the elements of the Redfield tensor \53 which are: 



R 



poo 

Jo 

k k ) 



(48) 



where 



4kiit) ^ mb\j){l\b^\k) + {^\b^\J){l\b\k)]^G^(t). (49) 



Here iG^iuj) = 7/(1 — e '^") is the Fourier transform 
of G>{t) and G<{lo) = G>(-a;), therefore, G<{t) = 
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G^{—t). This problem has a very interesting symmetry: 
the diagonal and off diagonal elements do not mix if the 
initial state of the charge particle is a diagonal density 
matrix in the BC. Therefore, we can divide the Redfield 
tensor elements in two parts, one which corresponds to 
population (diagonal elements) and other which corre- 
sponds to coherence (off diagonal elements). 

The Rii^nn elements which affect the population are: 



Here 



Ru,nn = / {x^niG> (r)e^'^"'^ + x^niG< (r)e-'""'^ } 
Jo 

= X^n[^G>iuJn^))]. (50) 

Here n ^ i and Xin — (|(n|6|i)p + |(n|6'l^|i) p), and 

Ru,u = ^^Xik[iG>{uj^k)] ■ (51) 

k 

Now wc can calculate the elements which are respon- 
sible for the coherence part. In the adiabatic regime we 
have 7~r2 — r2'<cri~ fl', i.e., where the BCs are not 
static and the mixing of pab and pcd (|47| as well as their 
conjugates cannot be neglected. Hence the non-zero el- 
ements of R tensor which affect the coherence are the 
following: 



7 



2„,,'\ 



2 

7 



Rcd,cd = [1 + + + lic^W ~ S'^W')] 

Rab.cd = — [1+5 -iw], 
Rcd,ab 



— [l-5-iw\ 



C : 

s - 

5- 

5' - 

w - 
I 

w -■ 



cos[(6'-6l')/2], 
sm[{e-e')/2], 

tea H" ^db-) 

Wca - Web, 
Wda - Web, 

^anhf^ 
2 V 2 



TT 



2tt 



and ip{z) is the digamma function. 



Now we can construct the explicit form of the matrix 
R which is: 



/ Ri.i 

2-0 







R2.1 









R3,i 

0' y- 







\Ra.i 



i?i,2 







zl 

R2.2 









R3.2 









i?4,2 



i?i3 

y+ 



OOP 



R2,3 



OOP 

OOP P 

POP P 

P R3,3 

P z+ 

P P 

POP P 

POP P 

P i?4,3 



RiA \ 



POP P 

POP P 

y; P 

i?2,4 

POO 

POO P 

POP P 

POP P 

OOP i?3,4 

OOP 

POO P 

POO P 

p z; p 

OOP i?4,4 / 



where z+ = [l + 0^6 + s^5' -f i{c^w - s^w')] , 

z- = -| [1 - c^5 - s^5' + i{c^w + s^w')] , 
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FIG. 5: Contour of the entanglement after time t = 5 for all 
values of Ej and e (Josephson junction Hamiltonian param- 
eters) The temperature in this case is equal to = 0, 7 = 1, 
K = vf-y = 0.45. 
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FIG. 6: Evolution of entanglement for the pure dephasing 
case, i.e, Ej = 0. We can see that whereas entanglement 
decay is exponential for the 1/f noise, it slows down for a 
bath of harmonic oscillator (inset) at zero temperature. 



y+ 



o 

c 7 



iw\ , 



— [l-5-iw\, 



«2JJ ■ 



Exponential of the matrix R is the matrix representa- 
tion of the evolution channel. Therefore, V = exp{Rt). 
This gives us the evolution for the qubit plus the charge 
impurity. From here we can find the evolution map act- 
ing on the qubit Vs (see appendix for more details). 

The two parameters e and Ej in the Hamiltonian for 
the charge Josephson qubit Hq play a crucial role in the 
decoherence properties of the system. For example, if 
Ej = 0, the system Hamiltonian Hq commutes with the 
interaction Hamiltonian. This situation is called non- 
demolition evolution or pure dephasing. In this case there 
is no energy exchange between system and the bath. On 
the other hand when we have e = 0, the system Hamilto- 
nian does not commute with the interaction Hamiltonan. 
Therefore, the two situations are qualitatively different. 
We present a plot of entanglement in the phase space of 
e and Ej by evolving a maximally entangled state of two 
qubits with the bath (of charge impurities) acting only 
on one qubit. The qubit is evolved for a fixed time t 
and entanglement calculated for different values of e and 
Ej, FIG. [SJ As seen from the figure, as the value of Ej 
increases, for fixed e, the entanglement in the system in- 
creases; on the other hand, for fixed Ej as e increases, 
entanglement decreases. This is counterintuitive because 
with the increase in the value of Ej, dissipation increases. 
FIG. ini compares the time-evolution of entanglement for 
the harmonic oscillator bath and the charge-impurities 
bath, both under pure dephasing. While entanglement 
decay is exponential for the case of 1/f noise, it is slower 
for a bath of harmonic oscillators. We compare, in FIG. 
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FIG. 7: Evolution of entanglement for 1/f (telegraph) noise 
at zero temperature. Here different curves represent the evo- 
lution of entanglement for different values of Ej (from to 1), 
with the curve at the bottom corresponding to Ej = and 
that at the top to Ej = 1, while e is fixed and equal to 1. 



[71 the time-evolution of entanglement for various values 
of Josephson energy {Ej ) starting with the pure dephas- 
ing case given by Ej =0. We see that with an increase in 
the value of Josephson energy, entanglement remaining 
in the system increases. This is consistent with FIG. [5) 
Decoherence produced by background charges depends 
qualitatively on the ratio k = v/j, where k <^ 1 denotes 
the weak-coupling regime and k > 1 is the strong cou- 
pling regime. The latter gives rise to qualitative new 
properties. We find (see FIG. |8(b")] ), for k > I, that the 
time-evolution of entanglement does not depend on n. 
This is in contrast to the weak coupling regime, where 
the time-evolution of entanglement does depend on k, 
as seen in FIG. 



8(a) where an increase in k leads to a 



decrease in entanglement. Naturally, decoherence due to 
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(a)Plot of the entanglement (concurrence) as a function of 
time (t) for different values of coupling strength k = v/'y, 
in the weak coupling regime, i.e, k <^ I. Here the range of 

K is from 0.05 to 0.5, with 0.05 corresponding to the 
uppermost curve, and 0.5 to the lowest (bottom) one. We 
can see that as we increase ft, entanglement decreases. 
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(b)Plot of the entanglement (concurrence) as a function of 
time (t) for different values of coupling strength k = v/j, 
in the strong coupling regime, i.e, ft > 1. Here the range 
of ft is from 5.05 to 5.5. We can see that in the strong 
coupling region all the curves converge. 

FIG. 8: Evolution of entanglement with respect to time, for 
different coupling strengths and temperature T — 0, Ej = 
1 = e. 



the bath forces entanglement to decay with time for both 
cases. 



1. Evolution operator with Bang-Bang interaction 

The Josephson charge qubit in contact with an 1// 
bath is now subject to fast pulses, under the Bang-Bang 
dynamical decoupling scheme. The Hamiltonian for this 
radio frequency pulse is the same as in Eq. (|40p . If the 
time for which a pulse is active is tt, then the evolution 
operator for the pulse may be written as = I ® iux 



where 2UTp = ±7r. The total evolution can therefore be 
written as 

Vtotal = {VpVs{T)f'^ (53) 

where 2Nt = t. Since the RF pulses act on the system 
for very short amounts of time, the evolution of the sys- 
tem can safely be assumed to be governed only by the 
dynamical map Vp for the time period during which the 
pulse is operating. As can be seen from the FIGS. IHlto 
[T2l with the application of bang-bang pulses, the system 
exhibits the phenomenon of ESD. 

Let us consider the case where Ej ^ e ^ 1. Let us 
also fix the pulse strength to he U = SOtt and ensure 
that the pulses act for very short times. As defined ear- 
lier, the ratio of the BC bias v and the switching rate 
7 defines the weak and strong coupling regimes, the for- 

V V 

mer designated by — ^ 1 and the latter by — > 1. 

7 7 
In FIG. [91 we plot the time-evolution of entanglement, 

with the coupling strength as parameter. For weak cou- 
pling, we find that Iesd initially increases with coupling 
strength. This continues till a turning point is reached 

V 

at — =0.38 when Iesd — 880. After this, with increase 

in coupling strength, tESD starts to decrease. As a re- 
sult, a kink appears in the corresponding entanglement 
vs time plot, FIGS.IHIini The receding of tESD with in- 
crease in coupling strength continues well into the strong 

V 

coupling regime, i.e. for 5.05 < — < 5.5. It, however, 

7 

does not go to zero, but rather chooses to saturate at the 

threshold value of tssD ^ 10, FIG. [iTl The "turning" 

and the "saturation" features are well captured in FIG. 
V 

I12[ where we plot tssD against — and keep the pulse 

7 

strength and durations fixed. We observe a crossover 

V 

phenomenon around — ~ 0.38, where the value of tESD 

rises sharply, only to fall back again even quicker. 

In FIG. [131 we plot the behavior oHesd with Ej and 
find that, as we increase Ej and thus move away from 
the pure dephasing situation, the time to ESD keeps in- 
creasing. 



VII. CONCLUSIONS 

The importance of the sustenance of entanglement in 
quantum systems cannot be overstated. In this paper, we 
have studied a variety of control procedures aimed at do- 
ing exactly that. Most of these are designed to suppress 
decoherence at the level of single qubits. A majority of 
the systems considered in this paper are qubits coupled 
with harmonic oscillator baths at finite temperature T, 
the couplings being either of dissipative or of the dephas- 
ing type. The time-evolution of entanglement when such 
a bath acts on one side of the two-qubit maximally en- 
tangled state is known In the commonly occuring 
dissipative case, entanglement decays asymptotically at 
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FIG. 9: Effect of bang-bang decoupling on entanglement, in 
the weak coupling regime. If we compare this plot with FIG. 
(8(a) I, with the curves corresponding to the same values of 



K, we see that the bang-bang decoupling causes entanglement 
to disappear faster in time for a fixed value of the coupling 
strength. Here the parameters are as in FIG. (8(a) land the 
pulse strength is U = 50n with time for which the pulse was 
activated is Tp = 0.01. In the inset we have the evolution of 
entanglement for very small range (0.01 to 0.1) of coupling k. 
The thickest curve is the one corresponding to At = 0.38. This 
curve is important in the sense that it has the largest Iesd- 




FIG. 11: Effect of bang-bang decoupling on entanglement in 
the strong coupling region. Here again we can see the effect 
of bang-bang decoupling on the entanglement if we compare 
this plot with FIG. (8(b) I. Here the parameters are as in FIG. 
(8(b) I and the pulse strength is f/ = SOtt with the pulse du- 
ration r = 0.01. 
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FIG. 12: Iesd is plotted as a function of coupling strength k. 
Here we can see that there is a clear distinction between the 
strong and the weak coupling region. As we increase n the 
tESD tends to freeze and asymptotic value of Iesd is around 
10. The parameters used are as in the previous plots. 



FIG. 10: Entanglement evolution for the coupling parameter 
range 0.3 < k < 0.4, with the uppermost curve corresponding 
to K = 0.3 and the lowest (bottom) curve corresponding to 
0.4. One can see from this plot the formation and disppear- 
ance of the kink. 

zero temperature, whereas it dies a sudden death at finite 
non-zero temperatures. Squeezing the initial bath states 
increases the time to ESD. 

As stated, the aim of most control procedures is to sup- 
press decoherence. For the case of photonic crystals, the 
design allows the system to conserve coherence when it is 
within the photonic band gap. Modulating the frequency 
of the system-bath coupling aims to suppress decoherence 



by shifting the system out of the spectral influence of the 
bath. In both these cases it is found that the suppres- 
sion of decoherence is accompanied by a corresponding 
increase in Iesd- 

However, it will be erroneous to naively suppose that 
this is the norm. The exact opposite phenomenon is ob- 
served for the case of resonance fluorescence, where the 
coupling between the bath and a two-level atomic system 
forced by an external resonant fleld, is modulated. It is 
seen that an increase in the external field frequency f2, 
the Rabi frequency, results in a faster decay of entangle- 
ment (FIG. [3]). A further non-trivial effect observed is the 
saturation in the time to ESD: Iesd does not go below 
a threshold value no matter what the Rabi frequency. 
In what could point towards a possible explanation of 
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Parcek for a useful discussion. 



Appendix (A): Calculation for concurrence 

Concurrence is a measure for entanglement of forma- 
tion for a mixed state of two-qubit system given by Hill 
et al. [1^. Concurrence is defined as C = max(0, Ai — 
A2 — A3 — A4), where {Ai} are the square root of the eigen- 
values of the matrix R = pp. Here p = 02 ® <J2P* 02 ® 02 
and the complex conjugate is taken in the standard basis. 

For some simple density matrices we can calculate the 
concurrence very easily. For example, consider the ma- 
trix 



FIG. 13; Iesd is plotted as a function of Ej. This shows that 
as we go away from pure dephasing the Iesd increases. The 
parameters used are as in the previous plots. 



this phenomenon, we observe that the sudden death time 
stops being dependent on the Rabi frequency at = — ; 
strikingly, this happens to be the boundary between the 
overdamped 17 < — and undcrdamped > — regimes. 

The latter and significant portion of the paper, com- 
prising Sec. VI, is entirely devoted to dynamic decou- 
pling schemes where RF pulses, applied at short time- 
intervals, smoothes out unwanted effects due to environ- 
mental interactions. We discuss two qualitatively differ- 
ent system-bath models: the first being the usual qubit 
and harmonic oscillator bath pair with pure dephasing or 
QND interaction; and the second being a bath of charge 
impurities, simulating 1// (telegraph) noise, acting on a 
Josephson-j unction charge qubit. In the former, the de- 
cay of entanglement is always asymptotic, the only effect 
of the RF pulses being in a slowing down of the rate of 
loss. For the later case, we see that though the entan- 
glement decay is exponential in the absence of quick RF 
pulses, an application of the pulses brings out the phe- 
nomenon of BSD. A very interesting phenomenon, ob- 
served in the strong coupling regime, was the decrease 
in the time to ESD with increasing pulse strengths. This 
is extremely counterintuitive, and brings into perspective 
the fact the in the non-Markovian strong coupling regime, 
the dynamics of entanglement can be different than that 
of decoherence. This calls for the need to have careful 
exhaustive studies of entanglement in these regimes. 
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The set {Ai} is equal to {(1 + e-T*)/2, (1 - e-''*)/2}. 
The concurrence for this state is C = e^^^ . 



Appendix (B): Calculation for the evolution map Va 
acting on the qubit from the evolution map of the 
qubit plus charge-impurity 

The matrix representation of the evolution map for 
qubit plus charge-impurity V is achieved by exp(i?t) 
where R is given in Eq. (|52|) . This map is in the basis 
{|^±)|*)} ® {|^±)b)}- To make it computationally eas- 
ier we need to write it in the basis {|0±)|6'±)} ® 
since the matrix representation of the map acting on the 
qubit is in the basis {|0±)|^±)}- To change the basis 
we need the assistance of a unitary matrix (in this case 
permutation matrix) P which is defined as: 
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P{\a)\b)\c)\d)) = \a)\c)\b)\d); 



|a)|6)|c)|d) = 



where 



ai 
a2 



ai 



a2 



\a)\c)\b)\d) - 

=4>P = I(8)p(8)I, 



/ 1 \ 

10 

10 

V 1 / 



(59) 



di 
d2 

di 
d2 



(60) 



(61) 
(62) 



After conjugating the matrix by P we get: 



V = PVP' 



(63) 



We can write this 16 x 16 matrix as a 4 x 4 matrix, 
where each of the element itself is a 4 x 4 matrix 
where «,j £ {1,2,3,4}. Then the map Vs acting on 
qubit is simply = Tr(Zy). From here we can get 
the corresponding M matrix. Solving it analytically, in 
the present case, is not easy. Therefore, we use numerical 
methods to calculate the evolution operator and entan- 
glement evolution for a system of qubits. 
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